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[57] ABSTRACT 

An ultrasonic wave is applied to an anisotropic sample 
material in an initial direction and the intensity of the 
ultrasonic wave is measured on an opposite surface of 
the sample material by two adjacent receiving points 
located in an array of receiving points. A ratio is deter- 
mined between the measured intensities of two adjacent 
receiving points, the ratio being indicative of an angle of 
flux deviation from the initial direction caused by an 
unknown applied load. This determined ratio is then 
compared to a plurality of ratios of a similarly tested, 
similar anisotropic reference material under a plurality 
of respective, known load conditions, whereby the load 
applied to the particular anisotropic sample material is 
determined. A related method is disclosed for determin- 
ing the fiber orientation from known loads and a deter- 
mined flux shift. 

1 Claim, 3 Drawing Sheets 
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METHOD OF DETERMINING LOAD IN 
ANISOTROPIC NON-CRYSTALLINE MATERIALS 
USING ENERGY FLUX DEVIATION 
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ORIGIN OF THE INVENTION 

The invention described herein was made by employ- 
ees of the United States Government and may be manu- 
factured and used by or for the Government for govern- 
mental purposes without the payment of any royalties 10 
thereon or therefor. 

This is a divisional of copending application(s) Ser. 

No. 07/797,507 filed on Nov. 22, 1991, now U.S. Pat. 

No. 5,209,123. 
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BACKGROUND OF THE INVENTION 

1. Technical Field of the Invention 

The present invention relates generally to load moni- 
toring and more particularly to a method of determin- 
ing loads and fiber orientations in anisotropic materials 20 
using energy flux deviation. 

2. Discussion of the Related Art 

In isotropic materials, the direction of the energy flux 
(energy per unit time per unit area) of an ultrasonic 
plane wave is always along the same direction as the 25 
normal to the wave front. In anisotropic materials, how- 
ever, this is true only along symmetry directions. Along 
other directions, the energy flux of the wave deviates 
from the intended direction of propagation. This phe- 
nomenon is known as energy flux deviation and is Ulus- 30 
trated in FIG. 1. An ultrasonic transducer U is coupled 
to an anisotropic crystalline material AM and directs an 
ultrasonic wave through the material as indicated by 
directional arrow N which is normal to the wave front. 
The anisotropic nature of the material causes the energy 35 
flux to deviate from the normal arrow N by an energy 
flux deviation angle A, resulting in an energy flux vec- 
tor V associated with the deviated wave flow W. The 
direction of the energy flux is dependent on the elastic 
coefficients of the material. This effect has been demon- 40 
strated in many anisotropic crystalline materials. In 
transparent quartz crystals, Schlieren photographs have 
been obtained which allow visualization of the ultra- 
sonic waves and the energy flux deviation. 

The energy flux deviation in graphite/epoxy (gr/ep) 45 
composite materials can be quite large because of their 
high anisotropy. The flux deviation angle has been cal- 
culated for unidirectional gr/ep composites as a func- 
tion of both fiber orientation and fiber volume content. 
Experimental measurements have also been made in 50 
unidirectional composites. It has been further demon- 
strated that changes in composite materials which alter 
the elastic properties such as moisture absorption by the 
matrix or fiber degradation can be detected nondestruc- 
tive^ by measurements of the energy flux shift. 55 

Graphite fiber-reinforced composites such as gra- 
phite/epoxy, graphite/magnesium and gra- 
phite/aluminum exhibit very high stiffness-to-weight 
and strength-to-weight ratios, making them excellent 
materials for lightweight aerospace structures. Since 60 
these structures are intended primarily to carry load, 
the in-situ and non-destructive determination of load 
and load induced quantities such as stress is very desir- 
able. The energy flux of stress waves, i.e., waves af- 
fected by a load, propagating through anisotropic crys- 65 
tals has been shown to deviate from the direction of the 
normal to the plane wave. However, there has been no 
known indication that this deviation can be correlated 


2 

with the amount of applied load. This lack of knowl- 
edge is not surprising since crystals are normally not 
used in load bearing applications. In addition, there is no 
known work concerning any energy flux deviation in 
highly anisotropic, non-crystalline materials such as 
graphite fiber-reinforced composites and consequently 
no known work concerning any effect of applied load 
on such deviations. 


OBJECTS 

It is accordingly an object of the present invention to 
determine load and load induced quantities such as 
stress and strain in anisotropic materials using energy 
flux deviations. 

It is another object of the present invention to accom- 
plish the foregoing object non-destructively. 

It is a further object of the present invention to ac- 
complish the foregoing objects in a simple manner. 

Additional objects and advantages of the present 
invention are apparent from the specification and draw- 
ings which follow. 

SUMMARY OF THE INVENTION 

The foregoing and additional objects are obtained by 
a method of determining a load applied to a particular 
anisotropic sample material according to the present 
invention. An ultrasonic wave is applied to the aniso- 
tropic sample material in an initial direction and an 
angle of flux deviation of the ultrasonic wave front is 
measured from this initial direction. This flux deviation 
angle is induced by the unknown applied load. The flux 
shift is determined between this flux deviation angle and 
a previously determined angle of flux deviation of an 
ultrasonic wave applied to a similar anisotropic refer- 
ence material under an initial known load condition. 
This determined flux shift is then compared to a plural- 
ity of flux shifts of a similarly tested, similar anisotropic 
reference material under a plurality of respective, 
known load conditions, whereby the load applied to the 
particular anisotropic sample material is determined. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a prior art device used to show the energy 
flux deviation which is dependent on the elastic coeffici- 
ents of the material; 

FIG. 2 is an apparatus according to the present inven- 
tion for determining loads in anisotropic non-crystalline 
materials using energy flux deviation; 

FIG. 3 is a representative graph correlating the signal 
ratio of two adjacent receiving transducers and the 
applied load; 

FIG. 4 is a graph of the energy flux deviation as a 
function of the fiber orientation angle in a no applied 
stress condition; 

FIG. 5 is a graph of the shift in energy flux deviation 
as a function of the fiber orientation angle due to a 1 
GPa load applied along the X3 axis, i.e., in the fiber 
direction; and 

FIG. 6 is a graph of the shift in energy flux deviation 
as a function of the fiber orientation angle due to a 0.1 
GPa load applied along the xi axis, i.e., in the laminate 
stacking direction. 

DETAILED DESCRIPTION 

Referring to FIG. 2, a transmitting ultrasonic trans- 
ducer 10 is acoustically coupled to a non-crystalline, 
highly anisotropic material AM, e.g., a graphite fiber- 
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reinforced composite such as graphite/epoxy, gra- 
phite/magnesium or graphite/aluminum. The ultra- 
sonic transducer 10 propagates an ultrasonic wave 
through the anisotropic material AM which initially has 
a wave front which is normal to directional arrow N. 5 
As discussed in greater detail below, an applied load LI 
and L 2 causes energy flux deviation angles of A°QL and 
A°QT. These propagation directions change with the 
respective applied loads. The directional arrows QL 
and QT are normal to the deflected wave fronts. 

The following variables are used throughout the de- 
scription of the invention: 

6 angle, measured from the normal vector N, of fiber 
orientation of anisotropic material AM; 15 

QT vector designating “quasi-transverse” particle 
displacement of the propagated ultrasonic wave; 

QL vector designating “quasi-longitudinal'’ particle 
displacement of the propagated ultrasonic wave; 

A°QT angle, measured from the normal vector N of * 
energy flux deviation in QT direction; and 

A°QL angle, measured from the normal vector N of 
energy flux deviation in QL direction. 

An array of receiving transducers respectively la- 25 
beled Vj, V 2 , V 3 , . . . Vjvare acoustically coupled to the 
non-crystalline, highly anisotropic material AM oppo- 
site transmitting transducer 10 . The receiving transduc- 
ers are arranged such that adjacent transducers can 
detect the desired increment of propagation shift, i.e., of 30 
the deviation of energy flux, of the wave. Such an array 
is disclosed in U.S. Pat. No. 4,499,770, the specification 
of which is hereby incorporated by reference. The 
transmitting and receiving transducers can be any con- 35 
ventional type such as piezoelectric or capacitive sen- 
sors. The receiver can alternatively be a scanned optical 
ultrasonic receiver or any array of such receivers. 

As known, the received ultrasonic wave is converted 
by each receiving transducer into a respective electrical 40 
signal having a voltage proportional to the wave inten- 
sity. In one embodiment, the shift in propagation caused 
by an applied load is monitored as a simple ratio be- 
tween two adjacent receiving transducers, e.g., V 2 :V 3 . 

A sample correlation between the applied load and this 45 
signal ratio is shown in FIG. 3. Material specific data 
curves are used to determine the applied load on that 
material with the measured signal ratio. 

In a more specific application, the effects of nonlinear 50 
elasticity on energy flux deviation in unidirectional 
composites such as graphite/epoxy (gr/ep) was consid- 
ered. The angle of the energy flux deviation, AQT or 
AQL, for such a material was shown to be a function of 
applied stress as a result of this elastic nonlinearity. This 55 
deviation was modeled as detailed below using acousto- 
elastic theory and second and third order elastic stiff- 
ness coefficients for T3 00/5 20 8 gr/ep previously deter- 
mined in NASA Contractor Report 4100, “Ultrasonic 
Characterization of the Nonlinear Elastic Properties of ^ 
Unidirectional Graphite/Epoxy Composites”, William 
H. Prosser, Oct., 1987. 

Assuming linear elasticity, the components of the 
energy flux vector (E ij) are a function of the linear elas- 55 
tic stiffness coefficients (Cyfc) and the spatial and time 
derivatives of the displacement vector (u/). The equa- 
tion is given by 


_ _ ( auk V aui ^ 

J - c v kl ^ ix m A it ) 

where the Einstein summation convention on repeated 
indices is assumed throughout this discussion. The angle 
of energy flux deviation can then be calculated as the 
angle between the energy flux vector and the normal to 
the plane wave front. 

To include non-linear elastic effects on the energy 
flux deviation, acoustoelastic theory is used. This the- 
ory predicts an “effective” linear elastic stiffness tensor 
(c that is a function of the second and third order 
elastic coefficients and the applied stress (cry), as de- 
scribed in a communication to applicants from Dr. 
David M. Barnett of Stanford University. The expres- 
sion for the “effective” stiffness tensor is given by 

c *nlij= kn!ij+ Vnfili ( 2 ) 

where 8 // is the Kronecker delta and k „nj is given by 


knlij — Cnlij + Crlijt-nr + Cnsifils + Cnlpfiip + Cnliq^jg + Cnlijuv^uv -f (3) 
Crlijuv^uv^nr + c nsijuv^uv^ls + c nipjuv£u)fiip + C n liq uv e uv €.jg. 


In this expression, Cjjkluv are the third order elastic stiff- 
ness coefficients and e z y are the strains resulting from the 
applied stresses. If the applied stresses are within the 
linear elastic regime, the strains are given by 

tipijkio-kl (4) 

where s ijki are the linear elastic compliances which are 
the inverse of the stiffnesses. 

Thus, if the linear elastic stiffnesses and compliances, 
the third order elastic stiffnesses, and the applied 
stresses are known, an “effective” stiffness tensor can be 
calculated. This can then be used to compute the 
changes in energy flux deviation as a function of applied 
stress which are a result of non-linear elastic effects. 

The effect of stress on the energy flux deviation was 
modeled for unidirectional T300/52Q8 gr/ep which was 
assumed to be transversely isotropic. The fiber axis was 
designated to be the X 3 axis while the laminate stacking 
direction which is perpendicular to the fibers was 
chosen to be the xj axis. The values of the previously 
measured non-zero independent linear elastic stiffness 
and compliance coefficients for this material are listed 
in Table 1. Likewise, the values of the previously mea- 
sured non-independent non-zero third order stiffness 
coefficient are given in Table 2. In both tables, the con- 
tracted subscript matrix notation is used. 


TABLE 1 


Linear Elastic Stiffness and Compliance Coefficients 

c u 

(GPa) 

Sij 

(GPa)- 1 

cn 

14.26 

sil 

0.092 

c 12 

6.78 

s 12 

-0.042 

C13 

6.5 

SI3 

—0.003 

c 33 

108.4 

S33 

0.0096 

C44 

5.27 

S44 

0.190 


TABLE 2 

Third Order (Non-Linear) Elastic Stiffness Coefficients 
c ijk (GPa) c ijfc (GPa) 

ci ii —196 ci55 —49.1 
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TABLE 2-continued 


Third Order (Non-Linear) 

Elastic Stiffness Coefficients 

c ijk 

(GPa) 

c ijk 

(GPa) 

C112 

• -89 

C 344 

-47 

CI 13 

-4 

c 133 

-236 

c 123 

65 

C 333 

-829 

C144 

-33.4 




Calculations were performed for elastic waves propa- 
gating in the xjX 3 plane. As in any anisotropic bulk 
material, three elastic waves will propagate along any 
direction in this plane. Of the three waves propagating 
in this plane, one of them is always a pure mode trans- 
verse (PT) wave with its particle displacement polar- 
ized perpendicular to the xi and X 3 axes, i.e., along the 1 
X 2 axis. The other two modes are quasi-mode waves 
with components of particle displacements both along 
their direction of propagation and perpendicular to it. 
One is a quasi-transverse (QT) mode wave while the 
other is a quasi-longitudinal (QL) mode wave. All three 20 
modes suffer energy flux deviation except for propaga- 
tion along the fiber axis (X3) and the laminate stacking 
axis (xi). These are symmetry directions along which all 
three modes are pure mode waves and none suffers 
energy flux deviation. 25 

The energy flux deviation was first computed for the 
condition of no applied stress. Then, calculations were 
performed for two different states of uniaxial stress. The 
first was stress along the fiber direction (X3) of a magni- 
tude of 1 GPa, as shown in FIG. 5. The second was 
along the laminate stacking xi direction with a magni- 
tude of 0. 1 GPa, as shown in FIG. 6 . These values are 
near the reported ultimate strengths of this material 
along the respective directions. This allows an estimate 
of the maximum effects of stress on energy flux devia- 35 
tion. For both conditions, the change in the energy flux 
deviation angle from the condition of zero applied stress 
was computed over the range of propagation directions 
of 0 ° to 60° from the fiber axis at 2 ° intervals. 

The flux deviation angles of the three modes as a 40 
function of fiber orientation at zero stress are plotted in 
FIG. 4. The fiber orientation angle is the angle between 
the fibers and the normal N to the wave front or the 
intended direction of propagation. A positive flux devi- 
ation angle implies the energy deviates away from the 45 
fiber direction toward the xi axis while a negative devi- 
ation means that the energy deviates toward the fibers. 
Over this range of fiber orientation angles, the energy of 
the QL and PT mode waves deviates toward the fibers 
while that of the QT deviates away from the fibers. 50 

In FIG. 5, the change in the predicted energy flux 
deviation due to the application of stress along the fibers 
is plotted as a function of fiber orientation angle. The 
energy of the QT mode wave suffers the largest shift in 
flux deviation, reaching a maximum of 3° at a propaga- 55 
tion direction of approximately 20 ° with respect to the 
fiber direction. The energy of the PT mode wave 
changes by a smaller amount in the opposite direction 
while the QL mode wave suffers a negligible shift. 

The relative magnitudes of the flux deviation shifts of 60 
the different modes can be explained qualitatively by 
considering the ratios of the magnitude of the non-linear 
elastic coefficients to the linear coefficients. The pri- 
mary elastic coefficients affecting the propagation of 
the PT and QL modes are those dominated by the ma- 65 
trix properties. These coefficients are cn, C12, C44, cm, 
and ci 12 . The magnitude of the non-linear coefficients 
are over an order of magnitude larger than the linear 
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coefficients in this case. However, the ratios of the 
non-linear to linear coefficients which dominate the QL 
wave (C33, ci 33, and C333) are much smaller even though 
the magnitudes of the individual coefficients are larger. 
Therefore, the effect of non-linear elasticity on the en- 
ergy flux deviation should be much smaller for the QL 
mode wave. The previous measurements of the effect of 
matrix degradation on energy flux deviation also 
showed a larger change in the flux deviation of the QT 
mode wave with almost no change in the QL mode 
wave. 

The shift in energy flux deviation due to applied 
stress along the xi axis is shown in FIG. 6 . Again the QT 
mode wave suffers the largest change in flux deviation 
angle while the QL mode is almost unchanged. It is 
interesting to note that the direction of the change in 
energy flux is in the opposite direction from the case of 
applied stress along the fiber direction. 

These calculations demonstrate the effect of nonlin- 
ear elasticity on the energy flux deviation of ultrasonic 
waves in gr/ep composite materials. The modes indi- 
cate the angles of fiber orientation and wave modes that 
suffer the maximum shift in flux deviation for the cases 
of applied stress considered. Although the models pres- 
ented were for bulk waves propagating through a thick 
composite material, the same effect is expected for plane 
plate waves propagating in thin plates. The longer prop- 
agation paths possible along plates would make the 
effect more measurable and thus could improve the 
stress resolution. This effect is used to develop a new 
nondestructive method for monitoring stress in compos- 
ite materials or as a new method for measuring their 
nonlinearity. 

In practice, data for flux shifts are generated for a 
specific material having a specific fiber orientation at 
known loads. A family of curves for different loads 
could be generated with each curve representing a dif- 
ferent load condition as shown by way of example in 
FIGS. 5 and 6 where only one load condition is de- 
picted. This data and a subsequently determined flux 
shift allows one to approximate the value of an un- 
known load which caused this flux shift. Specifically, an 
ultrasonic wave is applied to a particular anisotropic 
sample material in an initial direction and an angle of 
flux deviation of the ultrasonic wave front measured 
from this initial direction. As discussed previously, this 
flux deviation angle is induced by the unknown applied 
load. Next, the flux shift is determined between this flux 
deviation angle and a previously determined angle of 
flux deviation of an ultrasonic wave applied to a similar, 
or the particular, anisotropic reference material under 
an initial known load condition such as a no load condi- 
tion. This determined flux shift is compared to a plural- 
ity of flux shifts of a similarly tested, similar anisotropic 
reference material, or the particular sample material, 
under a plurality of respective, known load conditions, 
whereby the load applied to a particular anisotropic 
sample material is determined. The “similar” aniso- 
tropic materials) are similar in that they possess similar 
fiber orientation, i.e., they exhibit flux shifts similar to 
those of the particular anisotropic sample material. 

The unknown quantity could likewise be fiber orien- 
tation of a known material, wherein a determined flux 
shift caused by a known load allows for approximation 
of the fiber orientation, as shown in FIGS. 5 and 6 . An 
ultrasonic wave is first applied in an initial direction to 
a particular anisotropic sample material under a known 
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load and having unknown fiber orientation. Next, an 
angle of flux deviation of the ultrasonic wave front from 
the initial direction is measured. A flux shift is then 
determined between this measured flux deviation angle 
and a previously determined angle of flux deviation of 5 
an ultrasonic wave applied to the particular anisotropic 
sample material under a known load such as a no load 
condition. This flux shift is then compared to a plurality 
of flux shifts of a plurality of anisotropic reference mate- 
rials having a plurality of respectively known fiber 
orientations and similarly tested under the same load 
conditions as the particular material, whereby the fiber 
orientation of the particular anisotropic material sample 
is determined. These reference materials are “similar” in 
that they are the same composite as the particular sam- 
ple but have various fiber orientations. 

Many modifications, improvements and substitutions 
will be apparent to the skilled artisan without departing 
from the spirit and scope of the present invention as 20 
described in the application and defined in the following 
claims. 

What is claimed is: 
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1. A method of deteimining an applied load of a par- 
ticular anisotropic sample material comprising the steps 
of: 

applying an ultrasonic wave to a surface of a particu- 
lar anisotropic sample material in an initial direc- 
tion; 

measuring the intensity of the ultrasonic wave via an 
array of at least two receiving points located on an 
opposite surface of the particular anisotropic sam- 
ple material; 

determining a ratio between the measured intensities 
of two adjacent receiving points, the ratio being 
indicative of an angle of flux deviation from the 
initial direction caused by the applied load; and 
15 comparing this determined ratio to a plurality of 
ratios of a similarly tested, similar anisotropic refer- 
ence material tested under a plurality of known 
load conditions, wherein the load applied to the 
particular anisotropic sample material is deter- 
mined by correlating said ratios to previously de- 
termined ratios due to known loads in said refer- 
ence material. 
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